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Abstract. Nonlinear gyrotropic medium is a medium, whose nat- 
ural optical activity depends on the intensity of the incident light 
wave. The Kuhn's model is used to study nonlinear gyrotropic 
medium with great success. The Kuhn's model presents itself a 
model of nonlinear coupled oscillators. This article is devoted to 
the study of the Kuhn's nonlinear model. In the first paragraph of 
the paper we study classical dynamics in case of weak as well as 
strong nonlinearity. In case of week nonlinearity we have obtained 
the analytical solutions, which are in good agreement with the 
numerical solutions. In case of strong nonlinearity we have deter- 
mined the values of those parameters for which chaos is formed in 
the system under study. The second paragraph of the paper refers 
to the question of the Kuhn's model integrability. It is shown, 
that at the certain values of the interaction potential this model 
is exactly integrable and under certain conditions it is reduced to 
so-called universal Hamiltonian. The third paragraph of the pa- 
per is devoted to quantum-mechanical consideration. It shows the 
possibility of stochastic absorption of external field energy by non- 
linear gyrotropic medium. The last forth paragraph of the paper 
is devoted to generalization of the Kuhn's model for infinite chain 
of interacting oscillators. 

Introduction 

In the week light flux, that propagates through medium, polarization 
of electromagnetic field does not depend on the intensity of light and 
is defined unambiguously by the polarization of the field incident on 
the interface of vacuum and the matter under study This relationship 
is changed radically in nonlinear optics, where the indices of refraction 
and absorption are the functions of the radiation intensity. This may 
lead us to polarization multi-stability and polarization chaos. When 
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the intensity in the "input" is changed adiabatically in the "output" 
polarization is not yet the simple function and consists of stable and 
unstable areas [1]. At certain values of medium and radiation parame- 
ters, when we have in the "input" stationary radiation, in the "output" 
polarization undergoes auto-oscillatory changes or polarization pseudo- 
chaotic changes in time may happen, which have continuous frequency 
spectrum [1]. Polarization instability, due to nonlinear gyrotropy, is a 
matter of great interest for nonlinear optics of biological objects and 
exiton spectroscopy in solid states. The dependence of natural optical 
activity of isotropic medium on the intensity of electromagnetic wave 
is the simplest example of manifestation of nonlinear gyrotropy [2]. 

Nonlinear gyrotropy manifests itself clearly in such substances that 
are characterized by strong non-locality of nonlinear response. The ra- 
tio of characteristic size of medium to the wavelength plays the role of 
the spatial dispersion representative that determines the scale of non- 
locality. Therefore [3], the effects of spatial dispersion of nonlinearity 
and polarization instability connected with them may be important 
first of all in the cholesterol liquid crystals near the phase point of 
transition into the isotropic state (5/ X ~ 1CT 1 ) , in chiral biological 
macromolecules (S/X ~ 1CT 1 ) and in semiconductors near the absorp- 
tion resonances of exiton and bi-exitons (S/X ~ 10 _1 ) [1]. The idea was 
developed to explain the polarization instability, according to which 
nonlinear gyrotropic medium may be described as the unity of nonlin- 
ear oscillators presented by the Kuhn's nonlinear model [1]. The author 
of the mentioned article could obtain the implreit material equations 
of gyrotropic medium, in which cubic nonlinearity of molecular oscil- 
lators is taken into account and which permits the indistinct solutions 
for medium response. First Kuhn introduced one of the popular oscil- 
lator models for optically active medium. Its detailed description may 
be seen in [4]. And we concentrate on some details. The Kuhn's mol- 
ecule is formed by two oscillators, which are connected to each other 
by elastic interaction and they oscillate along the mutually perpendic- 
ular directions. The oscillators are separated in space. The Kuhn's 
molecule has right-left asymmetry i.e. it is a chiral molecule. The 
response of medium consisting of the Kuhn's molecules is non-local 
agitation of one of the oscillators conducted to the second end of the 
molecule. The size of the molecule determines the characteristic scale 
of spatial dispersion. In spite of simplicity the Kuhn's model is quite 
productive. For example, the ensemble of the Kuhn's molecules ori- 
ented randomly represents the simplest model of isotropic gyrotropic 
medium. If we assume that the distance between the oscillators is zero, 
or if we rule out the bond between oscillators, the ensemble of Kuhn's 
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molecules randomly oriented represents the simplest model of isotropic 
non-gyrotropic medium. Generalization of the Kuhn's model on the 
problems of nonlinear polarization optics presents one of the impor- 
tant directions of macroscopic modeling. This kind of generalization 
requires taking into account the nonlinear elastic forces in molecules. 
Let us examine the variant of nonlinear oscillator model according to 
the article [4]. The Hamiltonian of the molecule consists of kinetic 
energy of both charges, potential energy of elastic deformation and en- 
ergy of oscillators interaction taking into account nonlinear members. 
We restrict ourselves to the inharmonious of individual oscillator and 
to the case of cubic nonlinearity bond between oscillators: 

U (ui, u 2 ) = ^o( M i + u l) + + u t) + ^ u i u l + K u i u l + u i u l 

here uj is a natural frequency of individual oscillator, a determine there 
interaction, a, b, c determines anharmonicity of oscillations. 




a 



o 
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Fig. 1 The diagram of the Kuhn coupled unlinear oscillator model 
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Function U is invariant to the substitution of u\ by «2 and if the 
equilibrium positions rf^ and rf^ coincide with each other , mole- 
cule has two planes of symmetry and the chirality does not exist. 
But if rf^ 7^ r*^ an d vectors Ui,U2 and rf^ — are noncollinear, 
the Kuhn's oscillators system loses symmetry with respect to any fi- 
nite transformation and acquires spiral structure (chirality), i.e. it 
shall have different interaction with light of right and left polariza- 
tion. In the simplest case the direction of the oscillations is orthogonal 
ui = {x, 0, 0}, «2 = {0, y, 0} and equilibrium positions have coordinates 
(0, ± a Q /2) (a Q - size of a molecule). 

Having described briefly the Kuhn's model we remark that our aim 
is not the calculation of the optical parameters of medium. For us the 
Kuhn's model is a matter of direct interest as nonlinear model, in which 
dynamical stochasticity may be formed. Our purpose is to determine 
the parameters of the problem for which: 

1) system is precisely integrable; 

2) chaos is formed in the system. 

In spite of this success some questions are still open. The point is 
that, the Kuhn's generalized model [4] developed in [1] presents nonlin- 
ear dynamical system the study of which (not only in the case of weak 
nonlinearity, when we can use the methods of perturbation theory, as 
was done in [1]), is a matter of interest itself. The Kuhn's nonlin- 
ear model realized in nonlinear optics may be used for experimental 
observation of quantum chaos. 

The subject of this paper is to study classical and quantum nonlinear 
dynamics in the limits of the Kuhn's generalized nonlinear model. The 
paper consists of the following paragraphs: 

The first paragraph is devoted to study of classical nonlinear dynam- 
ics in case of weak nonlinearity, when we can consider the set of equa- 
tions as quasi-linear. We have obtained, using Van-Der Pole generalized 
method [5, 6], the analytical solutions. In case of strong nonlinearity, 
when it is impossible to use analytical methods, numerical methods 
are used for this purpose. We have studied correlation functions of 
the solutions and we have determined the value of Kolmogorov-Sinai 
entropy [7,8]. Fractal dimension is determined Using Grassberger and 
Procaccia algorithm [10-12]. 

In the second paragraph we have studied conditions of the Kuhn's 
model integrability. The exact analytical solutions written in the form 
of elliptical functions [13] are given. We have studied the case, when 
amplitude of external field was modulated by slow changing radio wave 
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field. We have studied the possibility of formation of chaos in a one- 
dimensional nonautonomous system using Melnikov's criteria [14]. It 
is shown that in case of the polychromatic pumping overlap of the non- 
linear resonances [15] occur and dynamical stochasticity was formed. 

The third paragraph was devoted to quantum-mechanical study of 
our problem. In particular we have studied the case when amplitude 
of external field was modulated with slow variable field. 

The last forth paragraph of the paper is devoted to generalization of 
the Kuhn's model for infinit clain of interacting oscillators. 

1. The case of weak nonlinearity 

As we mentioned in introduction, the aim of this paper is to study 
the dynamics of the Kuhn's nonlinear oscillations influenced by the 
plane polarized light wave. Let us assume that the frequency of the 
external field is Q then according to [1], the equations that describe 
the Kuhn's coupled nonlinear oscillations have the following form: 

eE 

x + uj 2 x + ay + Tx + ax 3 + b(y 3 + 3x 2 y) + cxy 2 = — - cos(fit + k z a Q /2), 

m 

eE 

y+uj 2 y+ax+Ty+ay 3 +b(x 3 +3y 2 x)+cyx 2 = — - cos(f2t— k z a /2), (1) 

m 

where u is a natural frequency, a, b, c are constants connected to an- 
harmonic oscillations, m and e are mass and charge of electron, E XjV 
are the transverse components of the external electric field, T is a phe- 
nomenological parameter, that describes dissipation, k z is wave vector 
projection, a Q characteristic of spatial scale of the given problem. Let 
us consider that the following conditions are held: 

— « 1, ^% « 1. 

Introducing dimensionless variables, with the help of the following 
transformation 

x y 
x — > — , y — > — , t — Ut 
a n a n 



"o 



equation (1) yields: 



where 



2 2 

x + w x + vi? y = ~^ X) y)X > + w cos ^' 

to 2 . a . uj 2 E v 



f(x, y, x) = -x + A x 3 + B Q (y 3 + 3x 2 y) + C Q xy 2 
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g(x, y, y) = -y + A y 3 + B (x 3 + 3y 2 x) + C Q yx 2 , 
aa 2 R _ X r _ A ( eE x ^ ' 2 

f(x, y, x) and ^(x, y, y) are nonlinear functions of the coordinates. When 
the coefficient fi placed before these functions is equal to zero fi = 0, 
the set of equations (1) becomes linear. When \i is not equal to zero 
but small, given set of equations is about linear one and that is why 
it is called quasi-linear. There are some methods of solution of quasi- 
linear set of equations. For example asymptotic method, based on the 
hypothesis of existence such kind of solutions that are similar in form 
to the solutions of linear set of equations. There exists method that 
does not require existence of these kinds of solutions, but is based on 
finiteness of any real system's bandwidth (so-called filter hypothesis). 
We have used the method of slowly varying coefficients, which employs 
averaging approach. First this approach was used by Van-Der Pole for 
solving oscillation problems [5] . Let us rewrite the set of equations (2) 
in the generalized coordinates in the following form: 

qi + A 1 q 2 + B 1 q 2 +nlq 1 = fif(q 1 , q u q 2 , q 2 )+D 1 sint+E 1 cost, ^ 
q 2 + A 2 q\ + B 2 qi+nlq 2 = jj,g(qi, q u q 2 , q 2 )+D 2 sint + E 2 cost. 

A 1 = 0, #i = nl = ^, fi= -1, D l = 0, E x = q x = x, 

Let us look for solutions of the set (3) in the form of solutions of linear 
set of equations: 

q± — a sin(A;it + j3±) + b sin(k 2 t + j3 2 ) + d\ sin t + e\ cos t, 

q 2 = aiasin(/c 1 t + fli) + a 2 bsm(k 2 t + f3 2 ) + d 2 sint + e 2 cost, (4) 

where a, b, f3i, f3 2 are slowly varying in time functions, k±,k 2 - the nat- 
ural frequencies of linear uniform set of equations, which are solutions 
of the following equation: 

ak A - (nl + n\- A X B 2 - A 2 B x )k 2 + n\n\ - B X B 2 = 0, 

(a = 1 - A X A 2 ) (5) 

and the rest of the coefficients are determined with the help of the 
following equations: 

A 2 k\ - B 2 nl- kl A 2 k\ - B 2 n\- k\ 



_ _ _ _ v 2 

^•1 9 79 a i 9 r~>> ^2 



n\-k 2 A x k\-B x ' 2 n\-k\ A x k\- B^ 
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, _ D^nl - 1) - D 2 (B 1 -A,) D 2 (n\ - 1) - D 1 (B 2 - A 2 ) 

di ^ , d 2 : 

E 1 (n 2 2 -1)-E 2 (B 1 -A 1 ) E 2 (nl-1)-E 1 (B 2 -A 2 ) 

e l — A J e 2 — a ) 



A = a - (n\ + n 2 - A X B 2 - A 2 B 1 ) + n\n\ - B X B 2 . 
Substituting these expressions in (5) yields: 

fj=l- A X A 2 = 1, 

,4_2^ 2 2 X_a 2 _ 
fc ft 2 ^ + ft^ ft^~ u ' 



(6) 



ujI . /cjf w 4 a 2 



ft 2 Vft 4 ft 4 ft 4 ft 2 ft 2 ' 

, 2 _ , ^ _ ^1 _ 
i_ fi 2 + fi 2 ' 2 ~ft 2 

— -B 2 — B 2 

ai = ^7^ = 1 ' a2 = ^ = - 1 > 

di = — = 0, d 2 = — = 0, 



A = 1 



2c 2 , c 4 a 2 _ (ft 2 - ul) 



2_,.,2\2_ a 2 



ei = 



ft 2 ft 4 ft 4 ft 4 

w 2 (a£ - ft 2 - E*a) 



A ft 2 - cj 2 - a 2 

1) - c, 2 (£> 2 -<F)-n) 

e2 



ft 2 — cj 2 — a 2 



o 



i?* = — 

-Ex 

As we have noted above, a,b,Pi,p 2 are slowly varying functions of 
time. Additional condition that will be imposed on these functions is 
the following: the time derivatives of the generalized coordinates must 
have the same form as in the case of the constant a, b, Pi, f3 2 coefficients. 
Taking into account the above mentioned, after time differentiation we 
obtain: 

qi = aki cos(fcit + Pi) + bk 2 cos(k 2 t + j3 2 ) — t\ sint, 
q 2 = a±aki cos(kit + Pi) + a 2 bk 2 cos(/c 2 t + f3 2 ) — e 2 sint, 
qi — a sin(A;it + f3i) + a(k\ + Pi) cos(kit + Pi) +b sm(k 2 t + p 2 )+ 
+b(k 2 + /3 2 ) cos(k 2 t + p 2 ) - ei sint, 
q 2 = aai sin(A; 1 t + Pi) + aa 1 (A; 1 + Pi) cos(/cit + Pi)+ 
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+ba 2 sin(k 2 t + p 2 ) + ba 2 (k 2 + $ 2 ) cos(k 2 t + j3 2 ) - e 2 sint, 
from which we get the following condition: 

asin(kit + Pi) +a$i cos(kit+Pi) + bsin(k 2 t + p 2 ) +bf3 2 cos(k 2 t + (3 2 ) = 0, 
aai sin(fcit + Pi) + aa±$i cos(kit + Pi) + ba 2 sin(k 2 t + p 2 )+ 
+ba 2 /3 2 cos(k 2 t + (3 2 ) = 0. 

Time differentiating of the time derivatives of the generalized coor- 
dinates qi,q 2 , substituting the obtained equations in the set (3), and 
taking into consideration (6) yields: 

q\ = aki cos(/cit + f3i) — aki(ki + Pi) sm(kit + Pi)+ 

bk 2 cos(k 2 t + p 2 ) - bk 2 (k 2 + $ 2 ) sin(k 2 t + (3 2 ) - e x cos t, 
q 2 = ak\oi\ cos(kit + Pi) — akiai(ki + Pi) sm(kit + Pi)+ 
bk 2 a 2 cos(k 2 t + p 2 ) - bk 2 a 2 (k 2 + $ 2 ) sin(/c 2 t + p 2 ) - e 2 cost, (7) 
q\ + Biq 2 + n\qi = /if(qi, qi, q 2 , q 2 ) + E x cos t, 
q 2 + B 2 qi + n\qi = fJ,g(qi, <ji, q 2 , q 2 ) + E 2 cost, 
aki cos(kit + Pi) + bk 2 cos(k 2 t + p 2 ) — akiPi sin(/cit + Pi) — 
-bk 2 p 2 sin(£; 2 t + p 2 ) = fif(qi, q u q 2 , q 2 ), 
aotiki cos(kit + Pi) + bk 2 a 2 cos(k 2 t + p 2 ) — akiPiai sin(A;it + Pi) — 
-bk 2 p 2 a 2 sm{k 2 t + p 2 ) = /ig(qi, qi, q 2 , q 2 ). 
If we introduce new variables £ = kit + Pi, i] = k 2 t + p 2 , then the set 
of equations for a, b, afti, b$ 2 is: 

d sin £ + aPi cos £ + b sin r\ + bp 2 cos r\ = 0, 
aai sin £ + aaiPi cos £ + 60:2 sin V + ba 2 p 2 cos 77 = 0, 
aki cos£ + bk 2 cost] - akiPi sin£ - bk 2 $ 2 sin r] = /J,f(qi, qi,q 2 , q 2 ), 
aaiki cos ^+bk 2 a 2 cosrj—akiPiai sin t;—bk 2 p 2 a 2 sin 77 = fj,g(qi, qi, q 2 , q 2 ), 
The determinant of this set is: 

sin£ sin?? cos£ cost] 

ax sin £ a 2 smr] axCos£ oncost? 
A;iCOs£ k 2 cosi] — A;xsin£ —k 2 sinr] 
ai/ciCos£ a 2 k 2 cost] — ax/^sin^ —a 2 k 2 sin rj 



A 



4:hk 2 



Ai = 



sin 77 
a 2 sin rj 
{if k 2 cos T) 

fig a 2 k 2 cos rj 



cos £ cos 1] 

ax cos £ a 2 cosi] 
— fcxsin^ — k 2 smi] 

-ax^isin£ — a 2 k 2 smr] 

= -2k 2 cos{f + g), 



= 2k 2 cos £fi(f+g) = 



A 2 = 



sin£ cos£ cost] 

aisin£ aicos£ a 2 cos?7 

fcicos^ fif — fcisin£ — A^sinr/ 

«i/cicos^ fig — ai/cisin^ —a 2 k 2 smr] 

= -2k 2 cos(f - g), 

sin£ sinr? cos 77 

cti sin £ a 2 sin r? a 2 cos 

A;iCos£ k 2 cosi] /i/ — /c 2 sin7? 

ai/cicos^ a 2 A; 2 cosr7 /ig — a 2 A;2sin77 

= 2A; 2 sin £(/ + (?), 

sin£ sinr? cos£ 

«! sin £ a 2 sin 77 cti cos £ 

fciCos£ k 2 cosr] — &isin£ /x/ 

aifcicos£ a 2 A; 2 cosr] — aiA;isin£ fig 

= 2ki sin 77(5 - /), 



= -2ki cosr) fi(g-f) = 



2k 2 sm£fi(f + g) 



2ki smr)fi(f - g) 



and therefore 



da 
~dl ~ 

db 

dt ' 

dfr 
a—— 
dt 

b dp 2 _ = , 
dt ' 



Ai 
A 

A2 
A 

A3 
A 

A4 
A 



2h 
1 



(/+<?) 



2fe 



if -9) 



cos£, 



cos 77, 



2^ / + ^. 
1 (<?-/) 



2£> 



sin£, 



smr?. 



The obtained equations are the set of equations (7) written in an- 
other variables. Suppose that a, b, Pi, l3 2 the variables change slower in 
comparison with the oscillations taking place in the initial dynamic sys- 
tem. If we average the obtained equations in periods 2n /k\, 2n/k 2) 2tt, 
then we obtain so-called cut equations for a,b,/3i,f3 2 variables: 



da 
~dt ' 

db 
dt 
d(3i 
~dT 
dfc 
dt 



1 

2k[ 
1 

2k 2 
1 

2k\a 
1 

2kol) 



(Fi + Gi) 
(F 2 - G 2 ) 
(F 3 + G 3 ) 
(G 4 - F A ) 



(8) 
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where 

2tt 2tt 2tt 2tt 2tt 2tt 

Fl = ~b>\ J J f QOS ^ dr l dt i F 2 = ^J J J f^dvdtdt, 

o o o o o o 

2n 2n 2n 2vr 2vr 2vr 

Fs = 4^3 / / / f sin Z d Z dr l dt i F4 = / / / f&nridridtdt, 

o o o o o o 

2tt 2tt 2tt 2tt 2tt 2tt 

Gl = ^J J J 9 c °s£d£dr]dt, G 2 = J J J g cos r}dr}d£;dt, 

o o o o o o 

2tt 2tt 2n 2n 2vr 2vr 

Ga = 4^3 / / / d^^drjdt, G 4 = ^J J J g sin r]dr]d^dt, 

o o o o o o 

Taking into account the values of a coefficient, let us write coordi- 
nates, there derivatives, / and g functions: 

qi = a sin £ + b sin 77 + ei cos £, 

g 2 = a sin £ — b sin 77 + e 2 cos i, 
51 = aki cos £ + 6/^2 cos 77 — ei sin t, 
g 2 = afci cos £ — 6A; 2 cos 77 — e 2 sin t, 

r 

/(?i> ?2, 9i, 92) = ^ ( a ki cos ^+bk 2 cost/— ei sint)+A (a sin £+6 sin 77 +ei cos t) 3 + 

+5 (a sin £—6 sin 77+62 cost) 3 +3-B D (a sin £+6 sin 77 +ei cost) 2 (a sin£— &sin7/+e2 cost) + 
+C Q (a sin £ + 6 sin 77 + ei cos t) (a sin £ — 6 sin 77 + e 2 cost) 2 , 

r 

<?(<7i, 52, 92) = ^( a ^i cos£— 6A; 2 cos 77— e 2 sint)+A (a sin£— &sin?7+ei cost) 3 + 

+B Q (a sin £ +6 sin 77+ei cos t) 3 +3-B (a sin £—6 sin 7/+e 2 cos t) 2 (a sin £ +6 sin 77+ei cos t) + 

+C Q (a sin £ — 6 sin 77 + 62 cos t) (a sin £ + b sin 77 + ei cos t) 2 

Substituting these functions and integrating yields the following re- 
sults for F and G functions: 

_ aVk x _ bTk 2 _ aTh bTh 

Fl ~ ~TT' ~ IT' Gl ~~ ~TT' ~ ~^T' 

F 3 = ^[6A b 2 -2b 2 c+3a\A +AB +C )+6(A +B )ej+A(3B +C )e 1 e 2 ^6B +C )e% 

F i = h -[QAy-2a 2 c+3b 2 {A -AB +C )+Q{A -B )e\+A{W -C ^ 

G 3 = ^[6A 6 2 -26 2 c+3a 2 (A +45 +C , )+6(A +5 )e 2 +4(3fi +C , )e 1 e2+(65 +C , )e|, 
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G 4 = ^[-QA a 2 +2a 2 c-3b 2 (A -AB +C ) + (QB -2C )e 2 1 +A(-3B +C )e 1 e 2 +(QB -2C )e 2 2 ]. 

Substituting the values of these functions in the set of equations (8)? 
yields: 

da aT db bT 

It = ~ TT' It = ~7T' ( ' 

^ = ^[6A 6 2 -26 2 C , +3a 2 (A +4fi +C )+4(3fi +C , ) ei e 2 +(A +6fi +C' )( e 2+e2)], 
at 0K1 

^17 = ^[QA a 2 +2a 2 C -3b 2 (A -AB +C )-A(3B -C )e 1 e 2 +(QB -C -3A )(e 2 1 +e 2 2 )}. 
at 0K2 

Separation of variables and integration of (9) yields: 

a = e~&, b = e""*, (10) 

A = ^[(-^){QA b 2 -2b 2 C o +3a 2 (A o +4B o +C o )}+A(3B o +C o ) ei e 2 + 



+(A + QB + C )(e\ + e 2 2 )]t, 
& = ^[(-^){QA a 2 +2a 2 C -3b 2 (A -AB +C )}-A(3B -C )e 1 e 2 + 

oK 2 

+ (6B -C -3A )(e 2 1 + e 2 2 )]t. 

As we can see from the results obtained, in the time interval < 
t p it is possible to neglect relaxation processes, and effects of non- 
linearity influence only the phases. In the time intervals greater then 
the previous one t > p dissipation processes become significant. Due 
to decrease of the amplitude of the oscillation the influence of the non- 
linear terms becomes insignificant. As a result, the forced oscillations 
of the linear coupled oscillators are obtained. This result obtained by 
qualitative analysis is in good agreement with the numerical solutions 
(see Fig. 2, Fig. 3). 
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D 5 10 , 15 20 2a 30 



Fig. 2. The plot of dependence of the oscillating system displacement 
(the nonlinear coupled oscillators) on the displacement's rate obtained 
by numerical integration of equation (1) for the following values of the 
parameters: 



Ao,= 0.3, B Q = 0.25, C = 0.25. 
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Fig. 3 The plot of dependence of the oscillating system displacement 
on time (the nonlinear coupled oscillators) obtained by numerical inte- 
gration of equation (1) for the following values of the parameters: 



n 2 



10, 



0.5, 



a 



q n 2 

Ao,= 0.3, B = 0.25, C = 0.25. 



1.5, E y = 2E X 



In case of strong nonlinearity it is impossible to use the method of 
slowly varying amplitudes and for description of system's behavior it is 
necessary to use numerical methods. In case of strong nonlinearity dy- 
namical stochasticity may appear in the system. It is necessary to take 
stochastic description instead of dynamic one. For description of pro- 
cesses taken place in dynamical system the following statistical quan- 
tities become important: Kolmogorov entropy [7], fractal dimension 
of strange attractor [10-12]. These quantities, statistical by definition, 
are closely connected with the mechanical parameters of the system, 
specifically with local instability of phase trajectories. Our purpose is 
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to determine those values of parameters of the given problem for which 
chaos is formed in the system described by the set of equations (1). Fig. 
4 and Fig. 5 show these values obtained by numerical calculations. 




10 20 



Fig. 4. The plot of dependence of the oscillating system displacement 
on time obtained by numerical integration of the set of equations (1) 
for the following values of parameters: 
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0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 
CO 



Fig. 5. Fourier image of correlation function obtained by integration 
of the set of equations (1) for the following values of the parameters. 

g =«>■£ = 0.6, £ = 1.5,^ = 2*. 

A Q , = 3, B Q = 2.5, C Q = 2.5. 



In order to prove existence of chaos it is necessary to study correlation 
function of solutions. For the definition of Kolmogorov entropy we use 
Fourier transformation of correlation function: 

G x (t) = (x(t + T) X (t)), (11) 

where 

T 

(x(t + r)x(t)) = lim / (x(t + r)x(t))dt, 

o 

is averaging in time r c is correlation time, which is connected with Kol- 
mogorov entropy through expression h ~ -. As numerical calculation 
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shows correlation length is equal to r c = ^ = 0.83 x 1CT 14 sec. The 
results of numerical calculations is shown on Fig. 6. 

In order to calculate fractal dimension of the system's phase space 
we use Grassberger-Procaccia algorithm [10-12]. 

The sense of this algorithm is the following. Let's suppose we obtain 
state vectors ensemble {xi, i = 1,2, ,7V} by numerical solving of the 
set of equations, that corresponds to successive steps of integration 
of differential equation. Choosing small parameter e we can use our 
results for evaluation the following sum: 



1 N 

c{£) = Jl m oc n(n^T) E e{£ - \ Xi - x ^ (12) 



[ x < 

where 9 is a step function 6{x) — < \ x > ' ^ LCC0I ^ n S to Grassberger- 
Procaccia algorithm, if we know C (e) we can determine strange attrac- 
tor's fractal dimension with the help of the following formula: 



D = Um ^£). (13) 
log{e) 



Let us calculate C(e) for different values of e and present the results 
in log(e) and log{C{e)) coordinates. Supposed dependence C(e) is a 
power function e D . So the obtained plot must be a line with slope D. 
The results of the numerical calculations are shown in Fig. 6. 
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Log (E) 



Fig. 6. The results of the numerical calculations in logarithmic 
scale. The results are obtained with numerical integration of the set of 
equations (1), using the formulae (12), (13) for the following values of 
the parameters. 

uj 2 r a cu 2 

^ = lO,- = 0.5,- 2 =l,- 2 =l.5,E y = 2E x , 

A Q , — 3, B = 2.5, C Q = 2.5. 



This plot enables us to calculate fractal dimension: 

D _ Alog(C(e)) _ . ^ 
Mog{e) 

2. The condition of the Kuhn's model integrability 



Our next task is to determine the conditions of exact integrability 
of the set of equation (1). It is easy to see that in the absence of 
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dissipation and external field, the Hamiltonian of the set of equations 
(1) has the following form: 

•2 -2 2 

H = y + y+y (a: 2 +y 2 )+axy+^(x 4 +y 4 )+b{xy 3 +yx 3 )+^x 2 y 2 . (14) 
Let us introduce new variables in order to separate them: 

x = q± + q 2 

. y = Qi - Q2 

As a result (14) yields: 

4 4 

H = q\ + ql + aql + + 2bq 4 + y 1 ~ otq\ + ?>aq\q 2 2 - cq\q\+ 

+ ^-2bq 4 + C -f + q 2 u, 2 + q 2 2 u, 2 . (15) 

It is easy to see that when the following condition c = 3a is kept, 
our equation is integrable. When this condition is kept, the separated 
Hamiltonian is written in the following form: 

ifi = q 2 i + Bqj + Aq\ = E u H 2 = q 2 2 + B*q 2 + A*q 2 2 = E 2 , (16) 

where 

a + Ab + c. 2 
B = , A = a + u 2 ; 

a-Ab + c 2 
B = , A = —a + u . 

In consequence we have the following differential equation for the 
variable q. 

dqi 



The roots of the expression under root are: 



, , A ± v'.l 2 • 



2B 

Taking into consideration this our differential equation has the fol- 
lowing form: 

d qi \ f m f 2 A + \J A 1 + ABE[\ ( 2 A - VA 2 + ABE[ 



5 < 25 + q] 01 25 (17) 
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Our aim is to find solutions of equation (16) in the form of elliptical 
function [13, 16]. Because of this let us introduce notations: 

A + ^A 2 + 4BE 1 _ ^A 2 + 4BE 1 - A 
71 ~~ 2B ' 72 ~~ W ' 

Taking into account these notations equation (16) has the following 
form: 

1 dqi 




_ , ,, . . . .. (18) 

y/B^f 2 dt ' ' " - - ' 

It is readily seen that equation (17) can be brought to the following 
canonical form: 



< ^=\l^-k' 2 ql){l+k 2 q 2 ) ) 
where the coefficients k and k' must meet the condition: 

k' 2 = 1 - k 2 . 

The latter expression in our notation has the following form: 

72 7i 
2B (A + y/A 2 + ABE 1 ) - 2B 



y/A 2 + ABE 1 - A A + y/A 2 + ABE 1 

This equation imposes the following condition on the value of energy: 

E 1 = 2B + V45 2 + A 2 . 

Let us evaluate the coefficient placed in front of time derivative of 
coordinate. 

1 1 1 



/ R A+-y/A 2 +4B£i yj A 2 +4BEi-A y[E[ 
V IB 2B 

In consequence the equation has the form: 

^^ = xj(l + k 2 q 2 )(l-k' 2 q 2 ), 
here k 2 = 2B k' 2 = , 2B 

A+y/A 2 +ABE 1 y/A 2 +4BE 1 -A 

Analogous calculations for coordinate q 2 give: 
E 2 = q 2 2 + B*ql + A*qi 
dq 2 



(]i v E 2 - Bql - Aql 



1 dq 2 
fW> dt 



= yj(l + k* 2 q 2 )(l-k'* 2 q 2 ), 
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where k* 2 = , 2B ' k'* 2 = , 2B * . 

A* + ^J A* 2 +AB*E 2 V A* 2 +AB*E 2 -A* 

Thus we obtain the following differential equations for q 1 and q 2 
coordinates: 

1 dqi 
y/E[ dt 
1 dq 2 



l + k 2 q 2 )(l-k' 2 q 2 ), 
y/(l + k"q*)(l-k»*q*). 



\[E 2 dt 

Separation of the variables and integration in the first equation gives: 

t 

[ ^ = = [ y/Kdt = yfet. 

J J(l + k 2 q 2 )(l-k' 2 q 2 ) J 



y/(l + k 2 q 2 )(l-k' 2 q 2 ] 

In consequence we obtain equation that determines the function qi(t) 
in the implicit form. Reversing this expression yields final solutions in 
the form of elliptical function [14, 16]: 



qi (t) = sd(k, V-EiO- 

In its turn given Jacob elliptical function may be presented in the 
following form: 

sd(k,t) = sn(k,t)/dn(k,t), dn(k,t) = a/1 — k 2 sn 2 (k, t), 

where sn(k,t) is a Jacob sine function [14-16]. By analogy for the 
second coordinate we have: 

sn(k*, y/E^t) 



q 2 (t) = sd(k\ y/E 2 t) 
Thus we have the following solutions: 
qi {t) = sd(k, y/Eit) 



dn(k* , y/Evt)' 

sn(k, y/E~it) 
dn{k, y/E[t) ' 



*( t) = ^,^) = g£^g. m 

Based on these solutions let us evaluate whether the neglect of dissi- 
pation term is correct or not. Take into consideration that the period 
of the elliptical functions dn(k, \Z~Et) and sn(k, y/~Et) is 4K(k), where 
K{k) is full elliptical integral [13, 14]. Compare this period T = 4K(k) 
to characteristic scale of time for which dissipation processes are con- 
siderable. Taking into consideration equation (9) dissipation can be 
neglected if the following inequality is met: AK{k) ^> p. 

It can be seen readily that our problem is integrable if certain condi- 
tions are met. What will be consequences of taking into consideration 
external field and dissipate term, will chaos appear in the system? We 
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can make this question clear with the help of Melnikov's criterion [15, 
16]. This criterion can be used to find out existence or absence of ho- 
moclinic structure, which is one of the interesting objects of nonlinear 
dynamics. Homoclinic structure can exist in phase space of dissipate 
as well as conservative systems. Existence of homoclinic structure can 
be used as evidence of existence of chaotic dynamics in a system that 
in turn implies existence of continuum of trajectories of random be- 
havior. Melnikov's criterion reads: suppose we have one-dimensional 
integrable Hamiltonian H(x,p) = E, where x = X(t) and p = P(t) 
are solutions of undisturbed equations x = Hp, p = —Hx; f, g are 
periodic functions of time consisting of external field and dissipation, 
e is a small parameter (x = Hp + ef(x,p,t), p = —H g + eg(x,p,t)). 
If autonomy of the system is breaking chaos is formed in the system 
when the following function changes sign 



A(0) = / [H x (X,P)f(X,P,t) + H P (X,P)g(X,P,t)]dt, (20) 



where 8 is an initial time displacement. Thus using formula (1), (19), 
(20) the question of formation of chaos in the system is brought to 
study the behavior of the following function: 



+oo 



— oo 



+oo 




— oo 




— oo 




(21) 



— oo 



— oo 



The plot of the function A(6>) is given on Fig. 7. 
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Fig. 7. Dependence of A(9) function on 9 initial 
time displacement. 

As can be seen from Fig. 7 this function A(0) does not change sign. 
This fact confirms that when there is certain symmetry in the system, 
chaos is not formed despite of strong nonlinearity and nonintegrability. 
Moreover, if the following condition is met AK(k) ^> p and it is possible 
to neglect dissipation, system is integrable even if it is affected by 
external field. In this case our Hamiltonian has the following form: 

H(x, x, t) = H (x, x) + H NL (x) + eV(x, t), 

where 

H = ^(x 2 + Ax 2 ), H NL = Bx A ; 

1 9 

eV(x,t) = — xcos(fit), e = — ; 

eV(x, t) = ex cos(r2t), e 1 

x, x - coordinate and momentum of electron, uj - natural frequency, B - 
coefficient of nonlinear term, Q -modulation frequency of variable field 
amplitude. We can consider interaction of electron with external vari- 
able field as small perturbation. With the help of the canonical trans- 
formation of variables x = (2J/^/]4) 1//2 cos#, x = — {2I\f~A) 1 ^ 2 sin#., 
we come to new variables action- angle (1,9). In consequence we ob- 
tain: H a = I\[A. Considering that the resonance condition takes place 
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9 ~ Q, let us average the Hamiltonian with the fast phase 9 [17]. 



3 



2tt 



^ V v A, 

Here averaging is equal to averaging in time 2ix / \H\ = 2-n/VL. This is 
small in this study, and 

[21 [21 1 

eV(x,t) = e\ — - cos 6> cos(fit) = eW— --[cos(0 + tit) + cos(6» - Clt)] = 

V ./i V ^1 ^ 

= ey(J)[cos(0 + Vlt) + cos(0 - Vlt)}, 
H(I,(p,t) = H^ L + eV(I) cosy?, 

H? L = VAI + H NL , H nl = ^(I/VA) 2 B, V(I) = (I/2VA)^. 

Here we introduce a slowly varying phase ip = 9 — Qt, that changes 
slightly in time~ 2tc/H. Because 9 = VAt this corresponds to the 
resonance condition f2 = \J~A~. 

Averaging with the fast phase 9 corresponds to averaging in little 
time 2tt/VA. p> changes a little at this time and that is why: 



sV(x, t) = eV(I, p) = e \j — cos p> = eV(I) cos ip, 



H(x, p, t) = H(I, <p) = H? L + eV(I, <p), 
H^ L = y/AI + H NL . 
The equations of motion have the following form: 



tp = u>(I)+e&V(I,(p). 

The resonance condition is met for certain value I a that can be calcu- 
lated from the same resonance condition: tp = uj(I q ), where 

but resonance value of the action I D is calculated from the equality: 

VA-tt + 3B^ = 0. 

A 

Deviation from the resonance action is AI — I — I , AI <C I , only in 
this case p is small phase. 
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Taking into consideration the expression of V(I, (p), the equations of 
motion have the form: 

/ = eV(I) sin(/9, 
<p = + e 9V gp cos (p. 

Expanding in a power series of ^- <C 1, and taking into account 
resonance condition u(I ) = 0, the set of equations takes the form: 

Ai = eV(I ) siny9, 



= U > AI + £ <mn 



ai 



cos (p. 

I = Io 



r-NJ 



Let us evaluate approximations carried out. We have V(I ) ~ H 

Here we introduce dimensionless parameter of nonlinearity a = (ff)^j- 
Neglecting nonresonance terms can be made if the following condition 
is met for nonlinear resonance width Auj = \JeV(I )uj' <C \[A and 
relative width: 



Auj yfeVjJ elVAcu' . — 

— P= ~ 1=— = \ = y/ea < 1 

y/A VA V A 

Combining two inequalities yj\ <C 1 and yfea < 1 we have yfe <C 
yfa <C l/y/s or £ ^ a ^ V £ - Taking into consideration the letter 
inequality, let us evaluate the second term in phase equation: 

edV(I)/dI eV(I )/I eVAI /I e^AI f fl 



uj'AI ' uj'AI ' uj'AI uj'AI /I ' V a ^ ' 

Then the equations of motion have the form: 

AI = eV(I ) simp 



<p = uj'AI, < 22 > 

and taking into consideration the results of this expansion, the Hamil- 
tonian has the form: 

uj' 

H = — (AI) 2 + Vcos<p, (23) 

where V = eV(I ),uj' = ^\ I=Io . It can be readily seen that the given 
Hamiltonian, called universal Hamiltonian, resembles the Hamiltonian 
of pendulum, which has mass 1/uj' and momentum AI and which is 
put in gravitational field with acceleration g ~ V. The quantity 1/uj' 
is measure of system's inertia: the more 1/u/, the more difficult the 
system to be taken out from the state of resonance with the help of 
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pumping and vice versa. Taking into consideration the fact that in 
classical mechanics the problem of pendulum is exact solvable, our 
problem is brought to exact solvable problem. Coming out from this 
fact we can write solution of action directly: 

AJ+ = y/(E + sV(I ))/iu'dn(^(E + eV(I ))/u't; k) E > eV(I ) 

(24) 

(with period 2K(l/k)); 

A/_ = y/(E + eV(I ))/u'cn(y/(E + eV(I ))/u't; l/k), E < eV(I ) 

(25) 

(with period 4K(l/k)), where E is an energy, cn and dn are Jacoby 
elliptical functions (elliptical cosine and delta), K(k) is full elliptical 
integral of the second order. When E = eV(I Q ) or k — > 1 these two 
solutions are drawn near each other, and they have the form of instan- 
ton: 

A/+ -> AJ_ -> y/2eV{I ))/u'/ch{y/2eV{I ))/u't). 

3. Quantum-mechanical consideration 

In the final part of previous paragraph we were able to show, that the 
study of the Kuhn's molecule influenced by external field is equivalent 
to that of the universal Hamiltonian when the potential of action has 
certain symmetry. The aim of this paragraph is quantum-mechanical 
study of the problem. After writing the stationary Schrodinger equa- 
tion 

HV n = E n V n (26) 

for the universal Hamiltonian of the system Kuhn's molecule + pump- 
ing 

H = -- + V, V = lcos2cp (27) 
dip 2 

we come to the equation coinciding with the Mathieu's equation. 

+ (E n - V(l, = 0, (28) 

where E n — * 8E n /h 2 u' is dimensionless quatity, / is dimensionless am- 
plitude of pumping, u/ = ^jp- is the derivative of nonlinear oscillation 
frequency oo(I) with respect to the action I. Characteristic peculiar- 
ity of the Mathieu-Schrodinger equation is the specific dependence of 
eigenvalues E n (l) and eigenf unctions ty n (ip,l) on parameter / [18]. On 
the plane (E, I), on which spectral indices are given, this peculiarity is 
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manifested in the interchangeability of areas of degenerated and non- 
degenerated states [19], Fig. 3, Fig. 4. The boundaries between these 
areas are passing through branching points of energetic terms E n {l). 
Let us suppose that pumping amplitude is modulated with slow vary- 
ing electromagnetic field. We can take into consideration influence of 
modulation if we make substitution in equation (28): 



where AJ is amplitude of modulation in dimensionless units, v- fre- 
quency of modulation. We imply, that little change l(t) can encompass 
certain numbers of branching points N on the right and the left side 
of the separatrix line (see [19], Fig. 3, Fig. 4.). Not going into details 
we only note that on (E, I) plane different eigenfunctions correspond 
to different areas. For particulars see [19]. Due to modulation of the 
parameter l(t) the system goes from one area to another passing the 
points of branching. Thus in case of certain conditions: 

1) amplitude modulated electric field acts on the system; 

2) The coefficients of the potential of interaction satisfy the following 
requirements c = 3a. 

Then our problem is fully analogous with the problem studied in [20] . 
In the mentioned article, it was shown that, at the expense of indirect 
passage, energy of external field is absorbed and quantum chaos is 
formed. On the basis of aforesaid we can conclude, that experimental 
observation of quantum chaos is expected when amplitude modulated 
field acts on nonlinear gyrotropic medium. 

4. Generalization of the Kuhn's model. The chain of 



Our next task is to generalize the Kuhn's nonlinear model on the 
chain consisting iV oscillators. Our aim is to present nonlinear gy- 
rotropic medium as ensemble of the nonlinear oscillators. The question 
we are interested in is the following: how does initial perturbation prop- 
agate through such medium and whether system has soliton solutions 
or not? As an adiabatic approximation the Hamiltonian consisting of 
N oscillators has the form: 



l(t) — l + Al cos vt, 



(29) 



INTERACTING oscillators 



H = ^(^+y^n+^4)+^[a^n+i+^X+i+H 




n n 



We study the problem with fixed boundary conditions: 
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It is readily seen that the equations of motion of the Hamiltonian 
(30) have the form: 

x n = -[u 2 x n + ax 3 n + a(x n+1 + x n ^ 1 ) + 

+ c { X n-l X n + x n x n +l 

) + b(3x n _ 1 x 2 n +x 3 n _ 1 )+b(x 3 n+1 + 3x 2 n x n+1 )}. (31) 

In order to perform further calculations we introduce envelope func- 
tion ipi(t) = (—l) % Xi(t) [21, 22], that is a function of coordinate i. We 
imply that envelope function is smooth function of coordinate i. It 
enables us to proceed from discrete variable ipi(t) to continuous one 
x = li, where / is lattice period. Using Taylor expansion of continuous 
variable function in series we have: 

i>{x ±l)= t)j(x) + *(x)(±Z) + ^ xx (x)l 2 + 

+ ^ xxx (x)(±l) 3 + . . . (32) 

Substituting formula (31) into (32), neglecting / 4 and terms of higher 
power, we have: 

ij tt + (uj 2 -2a)ij + (7a-8b)ij 3 -[a + (6b-2c)ij 2 }ij xx -(6b-2c)ijij 2 x = 0. 

(33) 

With the aim to simplify our expression, let us suppose that the system 
has certain symmetry In particular, let us assume that the coefficients 
of non-harmonic terms satisfy the following requirement: a = b = 
3c. In order to solve equation (33) we use rotate wave approximation 
(RWA) [21,22]. In particular, let us find the solutions in the following 
form ip(x,t) = ip(x)cosut and neglect the terms with multiple of us 
arguments. As a result equation (33) yields: 

(cu 2 -2a- u?)i, - h^ 3 - aij xx = 0. (34) 

Let us show that equation (34) has exact analytical solutions that are 
periodic function of variable x. For this purpose we multiply equation 
(34) on ip x and integrate the result: 

u 2 + 2a-u%\ 2 36 4 2 

H +5-^ +^ x = Ci. (35) 

a J oa 

We are governed by article [22] to analyze further equation (35). Fur- 
ther we normalize ip x function on its maximum value ip m = ip max , f(x) = 
^j^-. Then we rewrite this equation in the form that coincides with en- 
ergy conservation law of point particle of unit mass put in external 
field: 

H= f 4 + U{f)i (36) 
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where 
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U(f) = -^ 2 m (l-f)(C 2 + f) (37) 



and 

Ci = ( I — - + i ) < 2 = r (38) 



^ 2 + 2a-^ 2 36 2 \ 2 d 



solutions of equations (36) depend on the values of constant C 2 . 

It has been shown in [22] that when C 2 > 0, functions if)(x) are oscil- 
lating functions of x which vary between minimum ifj m and maximum 
—ipm values. The case C2 = corresponds to the separatrix solution, 
which is represented by the single localized breather with ip(x) — > 
at x — > ±00. The third case (C 2 < 0) corresponds to solutions that 
changes between nonzero values ip ma x and ip m in- 

In order to satisfy boundary conditions (zero displacement in points 
x = and x — N + 1), it is necessary to take the first case (C2 > 0) 
[22], because only in this case boundary conditions can be satisfied. 
Spatial period of this oscillation is given by this expression [22] 

if 



nsj = A/4 - <39) 

Zero boundary condition at point x = N + 1 is met automatically, if 
half wavelength A/2 is the solution of the following equation: (A/2n) = 
N + 1, where n = 0, 1,2...N. Dissipation relation (39) determines 
spectrum of u frequency function of n and ip m . 

If we express f x from equation (36) and substitute it into (39), after 
passing to new variable sin a = f, equation (39) yields: 

tt/2 



o 

Here we make substitution: (A/4(iV+ l)/2n and use notation r 2 = C 2 . 
Using the formulae (36), (40) we obtain the following expression for 
spatial profile function: 

arcsui(ip/il) m ) 

/x 1 f da / 1 -t \ 

*(t)= 1 / 7 • 2 7 2U/2 - ^ 

36 12 J (sin ct + r 2 ) 1 /^ 

v rn 



16a 



As we mentioned above, the solution of the breather type corresponds 
to condition r 2 = 0. This condition determines actually frequency of 
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breather as a function of amplitude: 

"I = ^ - 2a - (42) 

If we take into account the condition r = in expression (42), we 
obtain: 

tt/2 

1 f da 



x(f) = —== / , 0<x<oo. 

V 16a rn 



16a ^™ arc sin/ 

Integrating this expression we obtain finally: 



/ 35 

*Pb(x) = ip m cos h' 1 (d —i/j^x). (43) 

Based on formulae (42), (43) we can conclude, that in non-gyrotropic 
medium excitations, whose frequency of propagation ujb is determined 
by formula (42), during propagation remains spatially localized. 
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